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Abstract 

We revisit the BPS spectrum of the supersymmetric CP(A— 1) two-dimensional 
model with Zjv-symmetric twisted masses mi (I = 0, 1, ...,N — 1). A related issue 
we address is that of the curves of marginal stability (CMS) in this theory. Previous 
analyses were incomplete. We close the gap by exploiting a number of consistency 
conditions. In particular, we amend the Dorey formula for the BPS spectrum. Our 
analysis is based on the exact Veneziano-Yankielowicz-type superpotential and on 
the strong-coupling spectrum of the theory found from the mirror representation at 
small masses, |m/| <C A. We show that at weak coupling the spectrum, with necessity, 
must include A — 1 BPS towers of states, instead of just one, as was thought before. 
Only one of the towers is seen in the quasiclassical limit. We find the corresponding 
CMS for these towers, and argue that in the large-A limit they become circles, filling 
out a band on the plane of a single mass parameter of the model at hand. Inside the 
CMS, N — 1 towers collapse into A stable states. 
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1 Introduction 



The two-dimensional M = (2,2) CP(N — 1) model with twisted masses exhibits 
a remarkable similarity to certain four-dimensional gauge theories. The underlying 
reason was revealed in [HE] where the supersymmetric CP(N — 1) model was demon- 
strated to emerge as a low-energy world-sheet theory on the non-Abelian strings, see 
[31 IH dl 12]. The coincidence between the BPS spectrum of the CP(iV — 1) model 
and that of M = 2 SQCD in a quark vacuum observed in [5] received a natural 
explanation. Namely, these two spectra, in fact, describe the same states but from 
two different perspectives: viewed from the bulk and from the world sheet. 

Superalgebra in the CP(iV — 1) model is centrally extended, with the central 
charge containing a topological term and (N — 1) Noether charges. The former 
includes a canonical contribution and an anomalous one, which is especially obvious 
in the limit of the vanishing twisted masses (e.g. [6]). The contributions of the 
Noether charges are proportional to the twisted masses. 

The masses of the BPS-saturated states reduce to the corresponding central 
charges. However, not every state is realized in the theory. The BPS spectrum 
depends on the value of the twisted masses mi (I = 0, 1, N — 1). If \mi\ ^> A, we 
are at weak coupling. In this regime the BPS spectrum contains states character- 
ized by the topological charge T = ±1, with infinite towers of heavier states with 
nonvanishing U(l) charges built on them. Not all states possible algebraically are 
dynamically realized as stable states in the spectrum. If |m/| <C A, we are at strong 
coupling. The BPS spectrum of elementary kinks shrinks down to N stable states. 

A general and detailed discussion of the BPS mass spectrum in the CP (AT — 1) 
model with the twisted masses was undertaken in [7J [5] on the basis of the Veneziano- 
Yankielowicz-type0 superpotential [IT] augmented by the twisted mass terms [TJ, 
in conjunction with a quasiclassical analysis. A concrete implementation of the 
spectrum of the dynamically realized BPS states on both sides of the CMS (as well 
as the CMS itself) in the CP(1) model was constructed in (6]. For N > 2 in the 
case of .Zjv-symmetric twisted masses the issue of the BPS spectrum and the CMS 
separating the strong and weak-coupling domains was addressed in [T2] and, later, 
in [13]. The latter two studies were based on the Dorey formula which, as we will 
show in this paper, is by far incomplete. The advantage of the i^-symmetric twisted 

1 The Veneziano-Yankielowicz superpotential in the CP(N — 1) models with no twisted masses 
was originally derived, in terms of twisted supcrficlds, in 8 , 9 , 10.. 
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masses compared to a generic set is rather obvious. Generally speaking, we have iV— 1 
twisted masses implying the parameter space of the complex dimension N—l (in this 
case it is more appropriate to speak of the decay walls rather than curves). Starting 
from N = 3 and higher, with generic masses, explicit analytic determination of 
the BPS spectra and CMS, even if possible, ceases to be instructive. On the other 
hand, with the ^TV-symmetric twisted masses, the model depends on a single complex 
mass parameter, everything becomes simple, and the would-be decay walls reduce 
to a set of CMS in the complex plane of mo, a single complex parameter defining all 
twisted masses. On the physical side the advantage of the iJjv-symmetric choice is 
also obvious: with this choice the twisted masses preserve the Zn subgroup of the 
Z2N symmetry of the massless model, which is a remnant of the axial R symmetry. 
Physics of the transition from weak to strong coupling becomes transparent. 

Our goal here is to work out a complete solution of this question. Technically the 
problem is due to the fact that the Veneziano-Yankielowicz (VY) superpotential is 
a multivalued function and, in calculating the BPS spectrum, one needs to carefully 
analyze its branches in order to select those which can be dynamically realized. We 
add two crucial ingredients missing in [5]. First, in the small- \mi\ limit the BPS 
spectrum is known [H] from the Hori-Vafa representation [15]. In other words, the 
mirror symmetry is instrumental in finding the BPS spectrum. At strong coupling 
only N — l stable states exist. These are the states that become massless at the 
Argyres-Douglas points [16], at which vacua collide. In our case, the -Zjv-symmetric 
masses, all vacua collide simultaneously. The vanishing mass requirement at this 
point, in conjunction with small mass formula, enables us to resolve the ambiguity 
in the VY superpotential and determine the spectrum of masses of stable BPS states 
both outside and inside the CMS in its entirety. We arrive at the conclusion that the 
spectrum found in [5] and consisting of a single BPS tower (for kinks interpolating 
between a given pair of vacua) is not capable of satisfying all the above requirements. 
We argue that the theory must instead have N — l towers, each of which is described 
by its own U(l) quantum number. We find this natural in view of the fact that the 
global SU(iV) symmetry is broken down to U(l) Ar ~ 1 by the twisted masses. For each 
of these U(l)'s there is a tower of states arising from quantization. Only one tower 
is distinctly seen in the quasiclassical treatment, however, and that makes it special. 
That is the tower described in [5]. All others are not resolved: in the quasiclassical 
limit which corresponds to large |m;| and large U(l) quantum numbers all N — 1 
towers fuse with each other. 
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With the BPS spectrum in hands we pass to the second stage: determination of 
the curves of marginal stability in the same CP (A — 1) model with ^Ar-symmetric 
twisted masses. For each of the BPS towers there must be a curve on which the 
relevant states decay. Altogether, we deal here with N—l curves of marginal stability. 
We find that the curve corresponding to the special "quasiclassical" tower always 
passes through the Argyres-Douglas point. For this reason we will call it primary. 
The primary curve is the innermost curve on the complex plane of mo, inside which 
only strong-coupling states are stable. Other, secondary, curves are larger in size and 
are (typically) near-perfect circles. When passing from the weak-coupling domain 
into the strong-coupling one, the towers of states decay on these curves, one by one. 
Finally, we consider the large- N limit of the theory and show that in this limit all 
CMS tend to round circles, with radii in the interval 

1 < |m|/A < e 2 , (e = 2.71828...) (1.1) 

in units of the low-energy scale A. In the limit of very large N the curves of marginal 
stability completely fill this interval, forming a round band. 

Organization of the paper is as follows. In Sec. [2] we review the gauged formula- 
tion of the model and present the VY superpotential with the twisted mass included. 
In Sec. Owe consider predictions for the BPS spectrum at small |mo| following from 
the mirror symmetry. Section H] is devoted to the Argyres-Douglas (AD) point. In 
Sec. [5] we obtain the complete BPS spectrum both at strong and weak coupling 
based on the VY superpotential and the conditions outlined above. Construction 
of iV — 1 curves of marginal stability is the subject of Sec. El Section [8] presents a 
summary of our results. 



2 Exact Superpotential 

In the gauged formulation, the bosonic part of the Lagrangian of the N = (2,2) super- 
symmetric CP(A — 1) theory can be obtained from four- dimensional SQED with an 
iV-plet chiral superfield, by a dimensional reduction to two dimensions [TU1 [7] , 

£ = ^ [\*% + \9A + + |VX| 2 



+ tD(\n t \ 2 - 2(3) + 2jj<7-^ |nf . (2.1 

V 2 
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in the strong-coupling limit e — > oo. Here eo is the gauge coupling while m 8 are 
the complex twisted mass parameters which can be introduced p2] via a background 
vector field in four dimensions. Moreover, 2/3 can be viewed as a coupling of the 
corresponding sigma model (sometimes r = 2/3 is used). 

We should mention that physically the mass parameters are not given by the 
masses mi themselves, but, rather, by their differences mi — m^ (or by mi — m, where 
m is the average mass). This is clear from Eq. ( 12. ip as one can shift all masses by 
any value via a redefinition of a. Thus, one can always impose the condition 



N-l 



J2 m i = °- ( 2 - 2 ) 



With arbitrary choice of masses we break the global SU(iV) invariance of the 
model down to U(l) Ar ~ 1 . We are interested in a special case when the masses are 
assumed to preserve a discrete Z N subgroup of the anomalous U(1)r symmetry. To 
this end they must sit equidistantly on a circle, 

mi = m -e 2nil/N , Z = 0,l,...,iV-l, (2.3) 

with a single complex parameter mo defining all masses. The condition (I2.2p is 
automatically met. We will see that in the CP(iV— 1) model, the physical dependence 
of the theory is, in fact, on tuq . 

The theory (12. ip classically has N vacua, which can be seen as solutions with all 
rii but one set to zero, 

m = (o, o), k = o,...,N-i, 

a = m k . (2.4) 

Note that we choose to number both the masses and the vacua from to N — 1. 

The chiral sector of this theory is described by an exact superpotential of the 
Veneziano-Yankielowicz type [11]. For a theory with twisted masses the Veneziano- 
Yankielowicz superpotential was derived in [HI [91 [10l [7] , and is obtained by integrating 
out the n l fields in Eq. (12. ip . 

We ff (a) = -ira + i- J> - mj ) jln ^jp- - lj . (2.5) 
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The hat over a indicates that it is actually a (twisted) superfield. In passing from 
flZU to (E3]) we rescaled a, 

" - 75- < 2 - 6 ' 

In Eq. (12. 5 p r is the complexified coupling, 

Q 

r = %r + — , with r = 2/3 . (2.7) 

27T 

The ultraviolet cut-off scale can be traded for the dynamical scale 

p = Ae 27Tr/N , (2.8) 

where r in the exponent on the right-hand side corresponds to the same normalization 
point /i. Then, the first term on the right-hand side of (12 .5p disappears while /i in the 
argument of logarithms is replaced by A. The CP(A^ — 1) model is asymptotically 
free. 

Now, Witten's formula [10] for the positions of vacua is 

7V-1 

J] (a - m,) = A". (2.9) 

1 = 

This is the equation for critical points of W c s(a). With the set of masses ( 12.31) Eq. 
(12.91) implies that the vacua of this theory are at 



^J\ N + m$ ■ e 2mp/N , p = 0, N-l. (2.10) 

Altogether, we have N vacua, in full accord with Witten's index. 

If | mo | is taken to be large, then it dominates over A in (I2.10p . and the vacua 
take their classical values (12 .4p . 

a p w m p , p = 0, N-l. (2.11) 

In the future, for determination of the BPS spectrum, we will need the values of the 
superpotential in the vacuum. It is not difficult to get 

W cff (a p ) = - ^{Na p + Y, mj \n^jp-}. (2.12) 



2 Without a loss of generality one can always assume A to be real and positive. We will stick to 
this convention. A possible phase of A is absorbed in tuq. 
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Various choices of the branches of the logarithms above correspond to distinct values 
of the Noether U(l) charges. 

In deriving expression (12.121) we used the Witten's relation (12. 9ft . With our nor- 
malization, the general formula for the mass of an elementary BPS state reduces to 
the difference of superpotentials in two neighboring vacua, 

m BPS = | Weff^p+l) - Weff(cTp) I ■ (2.13) 

Note that the elementary kinks are obtained if one interpolates between the neigh- 
boring vacua. The result is the same independently of which pair of neighbors we 
pick up. This is due to the Z N symmetry of the model. 

With the ^jv-symmetric masses (12. 3p . the theory at quantum level retains the 
Zn symmetry, which the masses do not break. This symmetry manifests itself in 
the invariance of the spectrum to the choice of vacua in (I2.13p . From now on we 
will choose the vacua o"o and <j\ as representatives and focus on the masses of kinks 
interpolating between the two, 

m BPS = IWeffOO - WeffM| (2.14) 

(in bulky expressions, we will omit the absolute value sign, keeping it in mind; 
certainly, m BPS is always a positive quantity). 

In the quasiclassical limit |Am| ^> A the leading contribution to the mass is 
given by the dominant logarithm in expression (I2.14p . 



W eff (a 1 ) W eS (a ) ~ ^Am-ln^ = r • Am , (2.15) 

Z7T A 

where Am = m\ — mo- 

It is well-known that the Veneziano-Yankielowicz potential, being a multibranch 
function, is too ambiguous. The degree of ambiguity of expression ( I2.14p is deter- 
mined by the logarithms in ( I2.12p . and can be symbolized by a linear combination 
of the masses rnj with arbitrary integer coefficients 

(integer) j • rnj . (2-16) 

If all these multiplicities were physical, one would have a set of Z N states in the 
spectrum, which is certainly not what is expected. Selection rules need to be formu- 
lated in order to restrict the set of the BPS states that actually exist. We postpone 
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the formulation of these rules until further, while for now do a simple mathematical 
trick which reduces the amount of ambiguity present in Eq. (12. 14ft . Our goal is to 
turn (12. 1411 being a function of all masses and two vacua into a function of a single 
parameter mo- 
Let us pull out a factor of e 2n '^ N from each term in W e fr(ci) which originally looks 

as, 

vM^i) = - ^{n*i + E m ^ ln ^x^}- (2 - 17) 

3 

Both terms in Eq. (12. 17ft do contain this factor. This move turns o\ into o"o, while 
shifts the numeration of masses in the sum. To keep the numeration of masses in 
the sum consonant with the logarithms, we also pull out e 2m ^ N from the argument 
of the logarithm. This constant addition vanishes when summed with rrij, while 
G\ inside the logarithm again turns into ctq. Effectively one arrives at 

Weff(tri) = e 2m/N W cS (a ) . (2.18) 

The question mark in ( 12. 18ft reminds us that this is not quite the whole story. 
We were not very careful with the phases of the logarithms, and could have easily 
missed (and actually did) some 2iri. This owes to the fact that In ab = In a + In b 
only modulo 2ni. But such an omission can just as well be merged into the gen- 
eral ambiguity ( 12. 161) of the expression (12.14j) . We can fuse this ambiguity into a 
single linear combination i N ■ m, where N is an arbitrary constant integer vector, 
N = (no, i ) with integer rij and rh = (m , m^-i)- Now we rewrite 

(J2HD as, 

m BPS = U (m ) + iN-m, (2.19) 

with an explicit function 

U (m ) = (2.20) 



i- (e-/- - l) {jV^/< + A" + Yl In V< + ^ 



m 3 



The branch of \J§ is fixed as follows: if x is a real positive number, (i) y/x is a real 
positive number; (ii) In x is a real number. Here rrij is meant to be a function of tuq 
as well, via ( 12.31) . 

The important statements are, 
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In Eq. (12.191) all ambiguities related to the logarithms are referred to N; 



• Now [/o( m o) is a fixed single- valued function of the complex parameter tuq in 
a certain region of the complex plane. 



Let us briefly comment on what is meant. The BPS spectrum exists everywhere on 
the complex plane of mo, and both expressions (I2.14p and (I2.19P must in principle 
describe it. But for that, they need to be made unambiguous at least in some region of 
the complex plane of m®. We claim that function [/ (mo) is single- valued in a domain 
of the complex plane wide enough to unambiguously determine the spectrum. The 
latter will be described in terms of the vector N. Exactly how wide the domain of 
the parameter m needs to be will be discussed in Section [51 It will also become 
clear why Eq. (1 2 . 1 9 [) is more directly related to the determination of the spectrum 
than Eq. 11213)) . 

In this way we will arrive at our master formula. We will be able to determine 
such an expression for the spectrum which will be consistent both at large and small 
|m|. In particular, we will find that the prescription obtained in [5], 



is valid only approximately, in the quasiclassical limit, which corresponds to large 
\m\ and large excitation number n. The result ( 12.211) was derived by quasiclassical 
quantization of the time rotation of kinks in the U(l) factors of the global group, 
and, therefore, must still be valid as an asymptotics. Below we argue, however, is 
that the description (I2.2ip is incomplete. 

3 Mirror Treatment 

Now we gradually pass to the discussion of what is known about the strong-coupling 
regime. There are N BPS kinks. They can be seen in the mirror representation 



• Vector N in Eq. (I2.19p has a direct relation to the spectrum. 



(-71, 71, 0, 0), 



n G Z , 



(2.21) 



mm, 




+ 




(3.1) 
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with 



As was shown in 
\rrij\ <§C A, 



j 

one can determine the masses of all N kinks near the origin, 
N 



2tt 



[e 2 ^ - l) A 



i [rrij 



m) 



(3.3) 



where m is the average mass, vanishing in the case. So, to the linear order in 
the mass parameter, one has N kinks with the masses given by a large A term, and 
the splittings determined by rrij themselves, 

N 



-f- 



,2ni/N 



) 



1 A 



j 



0, N-l. (3.4) 



4 Argyres— Douglas Point 

If all masses rrij sit on a circle the corresponding vacua will also be forced to sit on 
a circle, see (I2.10p . Therefore, only simultaneous collisions of all N vacua a p take 
place in the theory at hand. The Argyres-Douglas points [16] correspond to 



cr,. 







This occurs whenever rriQ 



A , or, in the m,Q plane, 



■in 



AD 



A e i7T/N ■ e 2nil/N 



I 



0. 



iV-1 



(4.1) 



(4.2) 



Of these, the most convenient for us will be the two AD points closest to the real 
positive axis, 



m 



AD(I) 



Ae in/N 



and 



The crucial observation about the AD point is that one of the N soliton states 
(13. 4p becomes massless at that location in the complex m plane. Briefly, if all the 
vacua merge at the AD point, then Eq. (12. 14ft tells one that one of the kinks becomes 
massless 

m BPS = Weff(ffi) - W cfr (a ) = 0. (4.4) 

Here we quote this as a qualitative statement and render it more precise in Section [5j 
For now it is almost trivial to note that, although both superpotential functions here 
are multivalued, there exists a certain branch on which the above difference vanishes. 



AD(II) 



Ae 



-iir/N 



(4.3) 
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5 BPS Spectrum 



We now turn to the discussion of the BPS spectrum in detail, with an emphasis on 
the strong-coupling domain. Surprisingly, the conclusions obtained in the strong- 
coupling sector will allow us to make implications for the weak-coupling sector as 
well. We first collect the results known and trustworthy about the CP(1) model [6], 
as the simplest and well-studied case, and then increase N. 



5.1 CP(1) Spectrum 

There are two kinks in the strong-coupling sector with quantum numbers 

(T, n) = (1, 0) and (T, n) = (1, 1). 

Here T bears the conventional meaning of the topological charge, the constant that 
multiplies W c fr(o"i) — W c fr(o"o), and n e Z is the Noether charge, connected with the 
angle collective coordinate of the unbroken U(l). 

The formula for the BPS mass in the CP(1) theory is well-known. We recover it 
from our master equations (I2.19P and (I2.20p . 

A/mg + A 2 + m 



CP(1) 



7T 



m\ + A 2 — mo In 



^Jm\ + A 2 - m , 



+ iN-fh. (5.1) 



The Argyres-Douglas points here are m^ D = ± iA, and correspond to the vanishing 
square root. Let us trace how different states become massless at these AD points. 
Say, choose the branch of the logarithm such that at the point mo = iA the 
logarithm equals in. Then the kink N = (1, 0) becomes massless at this point. 
Let us move from the point iA to — iA along a large circle, see Fig [TR It is easy to 
show that the argument of the logarithm in Eq. ( 15. lft will also follow a large circle, 
starting and terminating at —1, see Fig [Tb. However, this contour arrives to — 1 
from under the cut of the logarithm, and, therefore, the phase of the argument of the 
logarithm changes from in to —in. As a result, the kink masses now shift compared 
to m = + iA point, 

m BPs\~ = im + iN-fh. (5.2) 

One observes, that it is a different kink which becomes massless now (it actually is 
the (T, n) = (1, 1) kink, as we will see in a moment). It is not difficult to see that 
for this kink N = (0, 1). 
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(a) 



(b) 



Figure 1: (a) The large-radius contour in the mo plane, starting at the AD point iA 
and terminating at the point — iA. (b) The same contour shown in the plane of z - 
the argument of the logarithm in Eq. (15.1 ft . 

Let us show that Eq. (15. ip is in one-to-one correspondence with the result from 
the mirror theory (13. 4p . At small mo, the argument of the logarithm in Eq. (15. ip 
just turns into 1, while the first term in the bracket tends to A, 

'"bps ~ 2A + iN-m. (5.3) 

7T 

This does indeed agree with the spectrum (13.41) . if we define the two states to be 

N = (1, 0) and N = (0, 1). (5.4) 

These are the two states of the strong-coupling regime. How does it happen that 
at weak coupling one has the whole tower of states while at strong coupling there 
are only two? Similar to what occurs in the Seiberg-Witten theory, the states of 
the weak-coupling sector decay on the curves of marginal stability [T9J [201 12] • Only 
two states (15.41) survive when crossing into the strong-coupling region, and those are 
precisely the states which become massless at the Argyres-Douglas points. 

The kinks (15. 4p therefore must be part of the weak-coupling spectrum. We must 
be able to write a general formula for the latter, keeping in mind the quasiclassical 



12 



asymptotic ( 12.211) . Indeed, if one allows iV to be of the form 

JV = (-n + 1, n) , n e Z, (5.5) 

then the two kinks ( 15. 4 p correspond to n = and n = 1. On the other hand, at 
large n Eq. (15. 5p does reproduce the quasiclassical tower ( I2.2ip . We stress that it 
is Eq. (15.51) that describes the exact spectrum. The latter was obtained based on a 
meaningful and single- valued formula (I2.19p . In this illustrative example of CP(1), 
we stayed within one sheet of the logarithm. 

We need to mention, that for CP(1) (and only for CP(1)) our results are not 
incompatible with those found in [5] and described by Eq. (12.21 j) . 

( — n, n ) , n G Z . (5.6) 

The reason is that the extra unity found in ( 15 .5p can be included into the logarithm 
in Eq. (15. ip . That will alter the sign of the expression under the logarithm, after 
which Eq. ( 15. ip will match with the analogous expression quoted in j5] precisely. 
Already at CP (2) we will find that this coincidence is not valid. 

To summarize, the spectrum for the CP(1) theory is given by equations ( 15. ip and 
( 15. 5p . where n runs through all integer numbers in the weak-coupling region, while 
it is restricted to n = or n = 1 in the strong-coupling domain. 

5.2 General criteria 

We can now formulate the requirements for the spectrum of BPS states in the overall 
region of the complex mass parameter mo : 

• Quasiclassical limit — the spectrum at large m and large excitation number 
n must reproduce the semiclassical result (12.151) . (12. 2 1 p : 

m BPS ~ ^- (mi - m ) ■ In + i n • (mi - m ) ; (5.7) 

• Argyres-Douglas point — the only states that survive when crossing from weak- 
coupling into the strong-coupling region are those A^ states which become mass- 
less at the AD points; 

• Mirror spectrum — the latter A^ kinks must reflect the spectrum given by 
mirror formula (13 .4p in the small m limit. 
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Having formulated the criteria for spectrum selection from the variety of branches 
of the logarithms of the Veneziano-Yankielowicz superpotential, we can now proceed 
to higher N. Of these, CP(2) will be the first nontrivial example. 

5.3 Domain of m 

At first we address the question of the relevant domain of variation for the parameter 
m . Needless to say, both Eq. ( I2.14p and (I2.19P are multivalued in the complex plane 
of tuq. The analysis of the actual complex manifold of m with all its branch cuts 
is a separate problem. We do not need such an extended analysis for our purposes, 
however. We can limit ourselves to the domain 

-7T < Arg(m v ) < n. (5.8) 

The exploitation of the criteria formulated in Sec. 15.21 relies on the possibility of 
using the formula ( I2.19P (i) in the neighborhood of the Argyres-Douglas point; (ii) 
in the neighborhood of mo = 0, and (iii) in the region of large |mo| - such that 
there are no discontinuities (i.e. branch cuts) between these regions. 

The second goal is the calculation of the curves of marginal stability, which in 
principle does require one to move over the whole complex plane of rriQ. One can 
use the physical 2tt periodicity in the 9 angle to argue that the spectrum of the 
theory must be identical in the N sectors { e 2mk / N ... e 2m ( k + 1 )/ N j f mQ (although 
there might be a monodromy between the sectors). The spectrum and the curves of 
marginal stability will therefore repeat themselves in all these sectors, and thus it is 
enough to build them in a sector 2n/N wide in Arg mo- 

We choose the sector between the two Argyres-Douglas points 

m£ D = Ae ln/N and m£ D = Ae~ ln/N , (5.9) 

see Fig. [2j On the one hand the AD points are a useful reference, since the (smallest) 
curve of the marginal stability with necessity passes through the AD point (see 
Section E]). On the other hand, similar to what we did in CP(1), we will be able to 
identify the states that become massless at the (at least these two) AD points. To 
this end it will be useful to be able to move from one AD point to another. 

Third, and perhaps most important, our function (I2.20p is free of the branch cuts 
in the region (15. 9p . This is, perhaps, the most drastic difference from the function 
in the right-hand side of Eq. (I2.14p which does have branch cuts in this region. 
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One might want to move beyond the boundaries of the domain (15.91) . for example, 
in order to approach other AD points. Although, as we argued, we do not need that 
for our purposes, we will partly be able to address the latter aspect. 

5.4 Spectrum in CP(2) 

The case of CP(2) is sufficiently straightforward so that one can analyze function 
U (m Q ) in detail, 



U (m ) = - i- (e 2 ™/ 3 - l) J 3^/m 3 , + + ^ mj In 




(5.10) 

There are three Argyres-Douglas points in CP (2), of which we will be interested 
in just two, see Fig. |3j The third one, m^ = — 1 in the units of A, is 

separated by a pair of branch cuts of {To (mo) from each of our primary AD points. 
The branch cuts start at infinity (i.e. in the semiclassical regime) and continue into 
the strong-coupling area, terminating inside the unit circle (it is to some degree a 
matter of convention where to terminate the branch cuts). Since we always prefer to 
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Figure 3: The Argyres-Douglas points in CP(2) theory (in the units of A). Two points e ±m / N 
lie within our area of interest. The third one is separated from our area by branch cuts of the 
logarithms, shown with black solid lines. 



move along circles of a large radius in order to stay away from the strong-coupling 
area, these latter cuts prevent us from smoothly connecting to the third AD point. 
Nevertheless, we will still have our say about that point. 

In between the AD points m^ D ^ = e m l 3 and m^ D ^ = e~ m ^ 3 there are no 
branch cuts for (15.101) . and one can smoothly commute between them. On the other 
hand, it can be shown that the right-hand side of Eq. (12.141) does have branch cuts 
in our area of attention, and therefore, although not impossible, when moving from 
one AD point to the other extra caution needs to be exercised. 

The function ?7o( m o) contains three logarithms, and it is easy to calculate its 
value at the AD points. At m^ ® = e l7T ^ 3 it reduces to 

U (e m/3 ) = -im = -*e in/3 . (5.11) 

Therefore, a kink = (1, 0, 0) becomes massless at this point. On the other hand, 
it is straightforward to calculate the BPS mass in the small m limit — the whole 
function (15 . 1 0[) only gives a zero-order contribution in the mass, the only linear term 
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being given by i N ■ rh = i mo, 

m BPS = - A fe 2 ^ 3 - ll A + zm . (5.12) 

This is precisely one of the kinks described by the mirror formula (13.4j) . 

We can now smoothly slide from the AD point e i7r / 3 to the one e _Mr//3 , along a 
large-radius contour (i.e. we first radially reach a circle of a large radius, and then 
sketch an arc extending clockwise to Arg mo = — n/3, after which, finally radially 
returning to the unit circle distance; see Fig. H] where such a contour is shown for 
the case of general N). One can show by tracing the corresponding contours for 
y/rriQ + A 3 — rrij that neither of the logarithms of Eq. (15.101) steps over a branch 
cut. At the lower AD point, one has 

Uo(e~^ 3 ) = -inn = -ie i7r/3 , (5.13) 

with the same numerical value, the logarithms just got shuffled around. One has the 
kink N = (0, 1, 0) becoming massless at this location. Now similar to ( I5.12p . we 
obtain 

m BPS = - — (e 2 ™ /3 - l] A + fmi. (5.14) 

This precisely matches the second of the kinks predicted by Eq. (I3.4p . 

Now how about the third kink? Since we know that the right-hand side of equa- 
tion ( I2.19P is smooth in our area of attention, and we know that it must account for 
the whole spectrum, we should be able to accommodate for the third kink with a 
certain choice of N. Indeed, if one allows N = (0, 0, 1), then in the limit of small 
mass one has 

m BPS = - — (e 2 ™ /3 - l) A + *m 2 , (5.15) 

which returns us the third of the kinks of Eq. (13 .4p . 

Now, generally, the third AD point, m^ D = — 1, is fenced from our area by 
branch cuts. What it in particular means, is that one cannot recover the value of 
m BPS in that point by smoothly sliding m from one of our AD points, say, e t7V ^ 3 
over to the third one. It does not exclude this, however, from being performed with 
proper accounting for the branch cuts. Even though, we can still calculate the value 
of (I5.10P in the third AD point, arguing that whatever the branch cuts are, they 
are responsible for bringing Uo(niQ D ^ lu ^) to the resulting form. We know the answer 
anyway, since, it is again obtained from (I5.12p by exchanging places of the logarithms, 
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and thus numerically giving the same result. Only now this quantity, e tn ^ 3 is called 
differently — 1712, 

Z7 (— 1) = -im 2 = -ie in/3 . (5.16) 

Since we did not follow any contour to connect this result to the Eqs. f 1 5 . 1 2 j) and 
( 15. 13ft . we further solidify Eq. ( 15. 16ft with the following remark. In this latter equation 
( 15. 16ft . had it been not right, we could only be off by a branch of a logarithm of ( 15.10ft . 
One can fix this branch, by approaching the third AD point from the quasiclassical 
area, m — 0q = — oo, and argue that the result (15. 16ft is precisely what one 
would find. 

Independently of this, we notice that Eq. ( 15.16ft indeed is responsible for render- 
ing the kink ( 15 .15ft massless at m^ D ^ III ' ) . We therefore have ascertained the strong- 
coupling spectrum of CP (2), which consists of three kinks 



and (5.17) 



We immediately come to the following conclusion. The first two states are part of 
the BPS tower of the weak-coupling region 

N = (-n + 1, n, 0), (5.18) 

with, correspondingly, n = and n = 1. This corresponds quasiclassically to the 
set of states f 12 .21ft . The extra unity present in Eq. ( 15 . 18ft can be explained similarly 
to the one in Eq. f 15 .5ft . Note that this unity could have equally well been placed 
into the second position of ( 15.18ft . and yet the expression would constitute the same 
tower. 

There is no way, however, that the third state in Eq. (15. 1 7ft can be found from 
Eq. f 12 .21ft . The third state is an entity that makes the CP (2) model qualitatively 
different from CP(1), as it must be part of another tower of BPS states, 

TV = n , l). (5.19) 

This sequence of states is completely new and invisible in quasiclassics! More exactly, 
quasiclassically, the two towers blend together, as the difference between n and n — 1 
is negligible at high excitation numbers. Furthermore, the contribution of the unity 
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at the third position on the right-hand side of Eq. (15.191) is similarly inferior to both 
terms in Eq. (15.71) . 



m BPS ~ — Am • In — - — + in -(mi — mo) + im2, (5.20) 
2ti A 

and so is not seen quasiclassically either. The occurrence of the extra tower and the 
extra unity in (I5.19P is a quantum phenomenon, and is the result of the fact that 
(any) BPS spectrum formula must describe the whole spectrum in any smooth region 
of mo — in our case, the spectrum of states given by the mirror formula (13.41) . and 
the region shown in Fig. [21 correspondingly. 

In summary, we found in CP (2) two towers of BPS states, 

N(i) = (-ri(i) + 1, n {1) , 0), n ( i) G Z, 

N( 2 ) = ( - "(2), ™(2), 1), n {2 ) e Z. (5.21) 

These towers merge with each other at large n, and thus are not distinguishable in 
the quasiclassical limit. In the strong-coupling region, we find three states 

(1, 0, 0), 

N = (0, 1, 0) and (5.22) 
(0, 0, 1), 

which form a subset of the above BPS towers: with = and = 1, and 
iV(2) with rip) = 0, correspondingly. 

5.5 Spectrum in CP(AT - 1) 

The crucial role which was played by equation (" 12 . 1 9[) in the above discussion is a 
reflection of similarity of Eq. (I2.20j) and Eq. (13 .4p . Indeed, the structure of the 
latter two equations is identical - - both formulas possess a factor of e 2 " K% l N — 1. 
In Eq. (I2.20p this quantity multiplies a figure bracket which in the small m limit 
turns into A and in this way matches its counterpart in Eq. (13 .4p . The outcome of 
this is that, even if the branches of the general expression Eq. (I2.14p were fixed by 
some method, Eqs. (I2.19P and (I2.20p would still play a more prominent role than the 
former. One other confirmation of this will be given further, when we discuss the 
curves of the marginal stability in Section El 
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Our discussion of the CP(iV — 1) theory would be in vain with generalization of 
our previous treatment of CP (2). We start with two AD points within our reference 
area, see Fig. |2l Choose m^ ® = e i7r//7V first. In order to determine which kink 
becomes massless at that point, one needs the value of Uo(mQ D ^). 

We find 

t/ (mo D(I) ) = -tm = -ie m/N . (5.23) 

This equality can either be obtained by attentively looking at the logarithms of 
^o(^o) as a function, without having to calculate anything, or by actual explicit 
summing of all the terms in Uo(rriQ D ®). We instantly find from Eq. (I5.23P that at 



^AD(I) _ ^/N j g ^ Q Jjj n ]j 



N = (1,0, 0) (5.24) 

that becomes massless. Tending the expression (I2.20p to the limit of small mo, we 
see that this kink has the mass 



- (e 2ni/N - ll A + im (5.25) 
In v ) 



m BPs — 

Zix 

near the origin. This obviously is one of the states of the spectrum ( 13 .4p seen in the 
mirror representation. 

The point m^ ® can be smoothly connected with rri^ ^ = e~ l7T ^ N , see Fig. [H 
It can be shown that none of the arguments of the logarithms in the function [To (mo) 
passes through a branch cut, which is what we mean by smooth. To calculate 
Uo{rriQ D ( II - ) ) one, however, can use Eq. f)5.23p . The value of function Uo(m ) numer- 
ically is the same in all AD points, while at the point mo = e~ tn / N we just endow 
it with a different name, 

U (mv D{u) ) = -imy = -te in/N . (5.26) 

This means that the kink 

N = (0, 1, 0) (5.27) 

becomes massless at m^ D ^ n \ In the limit of small m , the mass of this kink takes 

1 f e 2*i/JV _ A A + imu (5.28) 
In v ) 



the form 



3 2tt 
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which, again, corresponds to one of the states described by the mirror formula ( 13.41) . 
The above two kinks are part of the same tower of BPS states 

N = (-n + 1, n, 0), (5.29) 

which exists in the weak-coupling region. Equations f)5.24[) and (I5.27P represent the 
states with n = and n = 1, respectively. 

All other AD points, although maybe disconnected from our region of interest by 
branch cuts, can still be easily dealt with, since we know the value of Uo(mo) for any 
AD point, 

£/ ( m AD) = -ie™/ N . (5.30) 
For the extra N — 2 AD points, one has, therefore, 

U (m£ D ) = -im k , with k = 2, N — 1. (5.31) 

It does not matter which AD point produces what mass on the right-hand side, the 
only important thing is that all other masses with k > 2 pop up on the right- 
hand side. As it can be easily seen, in terms of the vector N, this corresponds to the 
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kinks becoming massless at the respective AD points. The index k effectively shifts 
the unity in Eqs. (15. 24ft and (15. 27ft further to the right. The full set of the states 
which become massless then looks as, 



N = (0. 0. 1 0). (5-32) 
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Now in the limit of small masses m we precisely obtain the full spectrum H3.4[) 
predicted by the mirror representation, 

m BPS = - — [e 27Ti/N - lj A + im k , k = 0, N - 1 . (5.33) 

In the weak-coupling spectrum, the states (15.321) belong to towers. The first two 
states are part of one and the same tower (I5.29p . while of the rest of the states each 
belongs to its own one. This agrees with a generic expectation to have N — 1 towers, 
according to the breaking of the global SU(iV) by the masses, SU(iV) — > U(l) iV— 1 . 
The spectrum (I5.2ip of the CP (2) theory is then obviously extended for an arbitrary 
N, 





(-n (1) + 1, 
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ri(N-i), 


0. 


0, .. 
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(5.34) 



with all ri(k) integer numbers. Obviously, in the quasiclassical limit, these towers 
reproduce the asymptotics (I2.2ip . and therefore satisfy all three criteria which we 
listed in the end of Section IBTTl Note, that the total number of states in the physical 
weak-coupling spectrum is (N—l) Z, much less then the total multiplicity Z N allowed 
by 0232J. 

Equations (I5.34p and (I5.32p are our main results for the spectrum of the CP(A r — 1) 
model in the weak- and strong-coupling regimes, respectively. They exhibit a drastic 



22 



difference with what was thought of the spectrum of CP(N — 1) earlier [5], when it 
was argued that only one tower of the BPS states existed (for a given pair of vacua, 
of course). We stress that the emergence of the extra N — 2 towers is a pure quantum 
effect which could not be anticipated in the quasiclassical theory. All of N — 1 towers 
of states blend and become degenerate in the quasiclassical limit, making it hard to 
resolve them apart. 

6 Curves of Marginal Stability 

The fact that there are N — 1 towers of BPS states means that there should be N — 1 
curves on which those towers collapse. Looking at Eq. (I5.34p one can tell that the 
first of them is special, just by its appearance — the unity stands in one row with 
rim (it is a mere matter of convention whether to write this unity in the first or in 
the second position). We will find that this tower is also special for an objective 
reason — namely, its decay curve will necessarily pass through the AD point, while 
those of the other towers will not. For the same token, it will also be the innermost 
curve. That is, inside this decay curve, only strong-coupling states (I5.32p exist. 

We provide an important technical remark on the graphical illustrations in the 
following discussion. We will choose to draw curves of marginal stability in the 
plane of rather than in that of m . Because of 27r-periodicity in #-angle, a 
curve sketched in the mo plane repeats itself N times in each of 27r /N sectors of the 
argument of mo- This way, drawing a curve in the tuq plane is as informative. We 
call for attention however, that when drawing multiple curves, we will have to do 
special rescaling in the plane, for the sake of fitting multiple figures in one plot. 

Having the spectrum of the theory at hand, it does not cost an effort to write the 
equations for the curves of marginal stability. Each such equation needs to rephrase 
the condition that one of the towers (I5.34p completely decays. Let us write explicitly 
the expression for the mass of a BPS state as per the spectrum (I5.34p . 

m BPS = U (m ) + in^-K - m ) + im k , k = 1, N - 1 . (6.1) 

Here for sake of convenience we redefined the U(l) charge nm in f!5.34p with 

-> ri(i) + 1 . 

In terms of the expression for the mass, the "tower" is given by the term 

i n( fc ) • (mi - m ) . 
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For the tower to decay, the remainder in the right-hand side of Eq. (16.11) must be in 
phase with the latter term, or 



Re Uo{mo) + %mk = 0. (6.2) 
mi — mo 

This is the CMS equation. We obtain N — 1 curves here by letting k run from 1 to 
N — 1. 

Now, we can make a few assertions on CMS before starting drawing them: 

• The curves (16. 2 p either do not intersect or they completely overlap; 

• The primary curve corresponding to k = 1 (and only this curve) passes through 
the AD point. 

The first assertion is seen from Eq. (16. 2p directly. If the curves with k = p and 
k = q happen to intersect somewhere, then 

— 2 G K 6.3 

mi — mo 

at that place. But this ratio does not depend on the absolute value of m , nor on its 
phase, so it will remain real along both of the curves, which for this reason will have 
to completely match. We will find that this does happen all along. 

The second assertion, although obvious as well, deserves more attention. To see 
that the curve k = 1 passes through the AD point, we rewrite Eq. (16. 2p as 

Re Uo{mo) + im ° = 0, (6.4) 
nil — m 

where the substitution mi — > mo in the numerator obviously does not change the 
condition (in fact, Eq. (16. 4 p is exactly what we would have obtained from the spec- 
trum (15.341) if we had not done the shift of nm above). But we calculated the value 
of Uo(m ) at the AD point m ^ D( ' I ' ) in Eq. (I5.23p . which shows that the CMS condition 
is trivially met there. Then the curve has to pass through all AD points. 

In the rriQ plane there is just one such point, and, since a p vanishes in it, we can 
expand the above condition in o"o/m in its neighborhood. We have, 

Re V -4 = 0, a = ^. (6.5) 

^ rN + 1 m v 1 

r > 
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Figure 5: The cusp of the primary decay curve at the AD point in the plane. 



This function, up to a constant, is actually known as the so-called Hurwitz-Lerch 
transcendent, and is a special case of the hypergeometric function. For our purposes, 
however, we only need the leading order term of it, 



Re a 



N+l 







for a C 1 



(6.6) 



Solving this equation gives us the angle at which the k — 1 curve passes through 
the AD point, 



N + 2 7T 

N + 1 ' 2 ' 

0, (for CP (1) only) 

N + 2 7T 



(6.7) 



N+l 2 



where, as customary, we measure the angle from the real positive direction counter- 
clockwise. Equation (16. 7ft tells us that for N > 2 there are two opposite angles, 
and the decay curve has a cusp at the Argyres-Douglas point, see Fig. As N 
becomes larger, the opening angle of the cusp increases, and ultimately, when N is 
taken infinitely large the curve becomes smooth. 

CP(1) theory is special, as its curve has three angles +120°, —120° and 0° instead 
of two. The third angle corresponds to the extra flat part [—1, 0] of the curve sticking 
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Figure 6: The curve of marginal stability in CP(1) theory (wig plane), Ref. [6]. 

out of the AD point, see Fig. |6j Otherwise, the CP(1) case is the simplest, and so 
we start the illustrative part of our discussion with this theory. 

6.1 The decay curve in CP(1) 

There is just one curve of marginal stability in CP(1) and it is the primary one. It 
represents a sharp boundary between the areas of weak and strong coupling spectra. 
This curve and its features have been discussed in the literature [6j [12], so we just 
merely reproduce it with Eq. (16.21) . Figure [6] shows the curve in the plane of TOq. The 
whole graph is presented in units of A 2 . We restate the known facts that the curve 
has a cusp at the AD point, where also an extra part [—1, 0] of the curve connects. 
All three lines meet at the AD point at an angle 120° with respect to each other. 
The real interval [—1,0] is the analytical solution to the CMS condition. There is 
nothing like that for any other CP(iV — 1) theory, all other curves are just curves. 
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Primary curve in CP (2) 
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Figure 7: The primary decay curve in the CP(2) theory (mjj plane). 
6.2 CMS in CP(2) theory 

The CP (2) theory features two curves - - one primary and one secondary. The 
primary curve, as we established, has a cusp at the AD point with the opening angle 
135°, see Fig. 

The second curve is a circle-like loop with a radius of approximately 361 A 3 (as we 
will see, it is a circle to a very good accuracy). This is a very large circle, compared 
to the primary curve which is of the size A 3 . In order to plot them both on the 
same graph, we rescale the radial direction of tuq by taking the iV-th root from its 
absolute value, while leaving the phase of as is, 

|m c | = Jm^] 1 ^, Argm c = Argm^. (6-8) 

The subscript c stands for compressed. Figure [H] shows the two curves in the rescaled 
rriQ plane. This compression will appear useful for the large N case. Such a transfor- 
mation, certainly, distorts the cusp making it less expressed. But, Figure [7] assures 
us that it is there, and Figure [5] tells us that it is there for all CP (AT — 1). 

The radius of the external curve, 361 A 3 is found to be rather large. However, 
in the compressed m c plane (Fig. |SJ), and, equivalently, in the m plane, the curve 



27 



Decay curves in CP (2) 
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Figure 8: Both decay curves in the CP (2) theory (compressed m c plane, see text). 

has the size approximately 7.12 A. We will find later, that the maximal radius of all 
curves (in the large N limit) is e 2 A, which amounts to 7.39 A. Thus, already the 
secondary curve of CP (2) nearly saturates the maximum size. 

6.3 Larger-AT theories 

For the CP (3) theory, the curves are shown in Fig. [9j Again, we plot the curves in 
the plane where trq was radially compressed to m c . The two outer curves in CP (3) 
overlap, as a consequence of Eq. (16. 3p . 

m-2 — m% = m\ — mo ■ (6-9) 

The cusp in the primary curve is still present, although is flattened in Fig. [9] due to 
compression. 

For CP (4) the curves are shown in Fig. [TUJ The overlapping curves are shown 
with portioned lines. The radius of the outer most curve is 7.33A, which is very close 
to the upper limit! Figures [TT] and [El as an illustration, show the curves for the 
CP(6) and CP(9). 
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Decay curves in CP(3) 




Figure 9: Three decay curves in the CP(3) theory (compressed m c plane, see Eq. (|6.8 



Decay curves in CP(4) 
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Figure 10: The decay curves in the CP (4) theory (compressed m c plane, see Eq. ([6 
External radius is 7.33 A. 
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Decay curves in CP(6) 
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Figure 11: The decay curves in the CP (6) theory (compressed m c plane). Outer radius is 
7.37 A. 

Decay curves in CP(9) 
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Figure 12: The curves in CP(9) theory (compressed m c plane). Outer radius is 7.02 A. 
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Let us briefly review the main features of these plots. We observe that all decay 
curves look as nearly perfect circles. We do not set the goal here to rigorously prove 
that they are, but we are able to show that the secondary curves for a few first 
theories are circles to a good accuracy and have no reason to believe that this is not 
so for larger N. More detailed analysis of their shape lies outside the scope of this 
paper. 

We know that the primary curves are not circular because of the cusp, and if in 
our figures they look round, this is just an illusion caused by the compression. For 
example, for CP (9) with N as large as 10, the cusp opening angle would be 164°, 
which although close to 180°, would still be quite noticeable. It is true, however, 
that at larger N even the primary curves turn into circles (that is, in any plane). We 
will discuss this in the next subsection. 

However, the secondary curves are circular even for small N. Using the (A/m ) N 
expansion of the CMS condition (16.21) it is possible to prove the following statement. 
If it is known, that a particular curve passes at a large distance from the origin at 
least at one point, then such a curve must be a circle as perfect as (A/m ) N . For 
example, in CP(2), Fig. El the outer curve is a circle with the accuracy about 1/360. 
This remark allows us to see that the secondary curves are round for a few starting 
CP(iV — 1) theories, just by looking at the location where the curves cross the real 
axis. For the primary curves this statement obviously does not apply since they pass 
through the AD point at a unit distance from the origin. 

For larger N, the k = 2 curves come closer and closer to the primary one, and 
one might suspect that they start losing their shape. For example, the radius of the 
k = 2 curve in CP(14) theory is only ~ 14 in units of A 15 . A deviation from a 
circular figure at the level of 1/14 would be quite noticeable. However, a different 
effect takes over, which helps the decay curves stay "fit". 

6.4 Curves in the large- AT limit 

We now consider the question of the form of the decay curves in the limit of large 
N. It turns out that the analysis greatly simplifies. All the curves turn into circles 
of a certain radius. 

As usual, the primary curves are a separate topic. Consider the curve in the 
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neighborhood of the AD point, where it is described by Eq. ( 16.51) . 

n rN + 1 

Re V = 0, (6.10) 

^ rN + 1 v ' 



r>0 



and discard the unity in the denominator. The sum then reassembles into a loga- 
rithm, 

Re ^log(l - 4r) =0, (6.11) 

where we have replaced a with its definition. As N goes to infinity, the vacuum cr 
approaches the value mo, and we can replace the ratio do/mo in front of the logarithm 
with unity. The logarithm itself can be transformed into 

This equation is trivially solved if (mo/A.) N is a pure phase. The curve then has to 
be a circle everywhere, not only in the region of validity of expansion (I6.10p . We 
thus have shown that at large N all primary curves tend to a circle of unit radius 
(in units of A^). 

A similar idea, but in a different realization is used to reveal the shapes of the 
secondary curves. In the region of space where \(Jq\ > \mo\ (in particular, at large 
positive mo) one can expand the CMS condition in 1/a. This gives, 

( rA r_i) 9 2k -1 

r>0 N 

where k, again, is the number of the curve. At large N we can again drop the unity 
in the denominator on the left-hand side, which turns the sum into a logarithm. Also 
we expand the sine in the right-hand side to the leading order in 1/A. We have 

1 - ln|er /A| = cos 2 — - ix . (6.14) 

N 

Replacing <jq with m with an exponential accuracy, we arrive at the formula, 

ImJ = e 1 - coslJ ^ l7T , k = 1, N - 1, (6.15) 



in units of A. Even though this formula has been derived in the assumption of large 
N, it qualitatively gives a reasonable answer even for N as low as three! 
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Decay curves in CP(ll) 
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Figure 13: The curves in the CP (11) theory (compressed m c plane). 

The qualitative features, which are obeyed in all CP(iV — 1) theories at large 
N are as follows. The curves come in overlapping pairs, as given by the cosine in 
Eq. (|6.15p . The minimum radius of the CMS is one, and is saturated by the primary 
curve — this fact we already know. The smallest secondary curves correspond to 
k = 2 and k = N — 1. Their size depends on N. Finally, the maximum size is 

m™ ax = e 2 , (6.16) 

measured in units of A, and is reached by the k = j" 1 or the k = — + ^ - 1 curves, 
depending on parity of N. For odd N, the largest curve does not have an overlapping 
party. 

Interestingly enough, although these qualitative results were inferred from the 
large- N formula (16.151) . they are still valid for small N theories! In particular, e 2 
seems to be the absolute limit for the size of all of the curves (the deviation is that 
this limit is not attained at small N, but some curves do come very close as we saw 
above). Also, the pairing of overlapping curves described by the cosine appears to 
be correct for any N. 

We illustrate the large N limit formula (I6.15P in Fig. [13j The external thin circle 
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Decay curves in CP(IOO) 
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Figure 14: The curves in the CP(IOO) theory (compressed m c plane). 



envelopes the overall \mo\ = e 2 size of the figure. The circles of radii determined by 
Eq. ( 16.151) are shown with thin dashed lines, which perfectly overlay the numerical 
curves. In fact, the latter formula has a good agreement with CMS curves already 
for N = 8, while, as we mentioned, in overall it shows the right tendency already 
for N as low as three. 

As we increase N, the decay curves fill in the whole interval 

\m \ E 1 ... e 2 . (6.17) 

Formula (16.15P predicts how the curves lay into this interval. As a concluding illus- 
tration, Fig. [141 shows the curves of the CP(IOO) theory. 

7 Prospects for the Spectrum in Quasi-Classics 

We have used strong-coupling techniques, such as the mirror symmetry, in order to 
establish the form of the BPS spectrum, 

m BPS = U (m ) + % ra (fc ) • (mi - m ) + im k , k = 1, N-l. (7.1) 
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At weak coupling, it must also be possible to see the features of this result, via 
semiclassical methods. At large m, function Uo(mo) gives the canonical logarithmic 
contribution as shown in Eq. ( 15. 7ft . while the mass difference mi — mo constructs 
the towers. These two contributions are expected. A very non-trivial prediction of 
Eq. ( 17. ip is the last term, which distinguishes the different towers. If one re-writes 
all relevant equations in terms of the mass differences, then 

3 

since the average mass is zero in our case. This shows in fact that the last term in 
Eq. ( 17. ip represents a fractional U(l) charge. 

The phenomenon of the occurrence of a fractional charge is well-known [22| [6j [23], 
[21], and is due to the presence of fermions. In particular, for CP(1) one has a half- 
unit contribution to the central charge [B]. We observe this one half by re- writing 
the last term of Eq. (17. ip in terms of the mass difference 

irrii = —im = -iAm, (7.3) 

which acts as a half-integer addition to nm. Therefore in the CP(1) model there is an 
agreement between our spectrum and the known quasiclassical results. Recovery of 
the weak-coupling spectrum in a general CP(A — 1) theory via semiclassical methods 
is an important and non-trivial problem deserving a separate study. We plan to 
address this topic elsewhere. 

8 Conclusion 

We demonstrated that the weak-coupling spectrum of the CP (A — 1) theory with 
twisted 2^ masses is considerably richer than was thought before. In particular, 
there are A — 1 infinite towers of the BPS states. Our analysis relied on three 
important facts known about CP (A — 1): the strong-coupling spectrum only includes 
A states; these states become massless at the Argyres-Douglas points; and the 
quasiclassical spectrum contains a tower of states with integer U{1) charges. The 
knowledge of the exact superpotential is not sufficient as it is ambiguous. However, 
if one fixes this ambiguity at least in some region of the parameter space, then the 
superpotential must still describe the states in the strong-coupling domain. It is 
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possible to trace these states all the way from the weak-coupling region through the 
AD points and into the strong-coupling domain, where their masses can be fixed by 
comparing with the mirror representation. We emphasize that we observe N — 1 
towers instead of just one, which is naturally explained by the fact that the global 
SU(iV) symmetry is broken down to N — 1 copies of U(l), and the central charges 
of the model include terms proportional to N — 1 Noether charges. However, only 
one of these towers is seen quasiclassically. Furthermore, these N — 1 towers blend 
together in the quasiclassical limit, making it hard to anticipate their existence from 
the semiclassical analysis alone. 

Having obtained the spectrum, it is easy to construct the curves of the marginal 
stability. We find N — 1 such curves, one per each BPS tower. One of the curves, 
which we refer to as primary, is special as it passes through the AD point, with 
necessity. Inside this curve, only N stable states of the strong-coupling spectrum 
survive. We also considered the large- iV limit and argued that all curves tend to a 
circular shape, with the radius lying in between 1 and e 2 , in units of A. 

We note that we only analyzed the decay curves of elementary states. In principle, 
nonelementary kinks have their own series of curves [25]. Also, the theory must have 
the fermion-soliton bound states [26], for which there will exist CMS as well. 

As was shown in [5], the BPS spectrum of dyons (at the singular point on the 
Coulomb branch in which N quarks become massless) of M = 2 four- dimensional 
supersymmetric QCD with the U(iV) gauge group and Nf = N quark flavors, identi- 
cally coincides with the BPS spectrum in the two-dimensional CP(A— 1) model. The 
reason for this coincidence was revealed in [H [2]. The confined 't Hooft-Polyakov 
monopoles of the four- dimensional theory are represented by junctions of two differ- 
ent non-Abelian strings. The effective theory on the world sheet of the non-Abelian 
string is the two-dimensional CP(iV — 1) model [31 IU [U 12] • Confined monopoles of 
the bulk theory are seen in the world-sheet theory as kinks of the CP (A — 1) model. 
This ensures the coincidence of the BPS spectra in both theories, see [TJ [2] for more 
details. 

The above coincidence implies that we can use the results for BPS spectrum of 
the CP(iV — 1) model obtained in this paper to construct the physical spectrum 
of the BPS dyons in M = 2 supersymmetric QCD in the particular vacuum in 
which N quarks become massless. More specifically, the masses of the BPS states 
in four dimensions are given by the periods of the Seiberg-Witten differential, and 
have exactly the same form as the central charge of the sigma model (12.1 3D written 
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in terms of the Veneziano-Yankielowicz superpotential. The masses of dyons are 
determined by the contours that encircle the branch points of the Seiberg-Witten 
curve. The spectrum that we derived in the sigma model gives the prescription on 
how to build contours that correspond to stable BPS states. In particular, our result 
(I5.34p for the weak-coupling spectrum of the CP(iV — 1) model suggests that the 
dyon electric charges in the bulk theory, besides the towers built on the roots of 
the gauge group also contain contributions determined by its weights, cf. [27]. The 
question of the correspondence between the spectra of the BPS states and the decay 
curves of dyons of the two theories calls for a special investigation. 
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